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ABSTRACT 

We perform an extensive suite of N-body simulations of the matter power spectrum, in- 
corporating massive neutrinos in the range = 0.15 — 0.6 eV, probing the non-linear 
regime at scales k < 10/iMpc~^ at z < 3. We extend the widely used HALOFIT ap- 
proximation to account for the effect of massive neutrinos on the power spectrum. 
In the strongly non-linear regime HALOFIT systematically over-predicts the suppres- 
sion due to the free-streaming of the neutrinos. The maximal discrepancy occurs at 
k ~ l/iMpc~^, and is at the level of 10% of the total suppression. Most published 
constraints on neutrino masses based on HALOFIT are not affected, as they rely on 
data probing the matter power spectrum in the linear or mildly non-linear regime. 
However, predictions for future galaxy, Lyman-a forest and weak lensing surveys ex- 
tending to more non-linear scales will benefit from the improved approximation to the 
non-linear matter power spectrum we provide. Our approximation reproduces the in- 
duced neutrino suppression over the targeted scales and redshifts significantly better. 
We test its robustness with regard to changing cosmological parameters and a variety 
of modelling effects. 

Key words: neutrinos - cosmology: large-scale structure of Universe - cosmology: 
dark matter 



1 INTRODUCTION 

A variety of solar, atmospheric, reactor and accelerator neu- 
trino experiments have firmly established the existence of 
neutrino flavour oscillations, implying that neutrinos have 
at least two non-zero mass eigenstates ( Becker-Szendy & et 
al.|1992[[Fukuda fc et a l. 1998; Ahmed fc et al.|2004[ ). These 
experiments yield a lower limit on the sum of the neutrino 
masses; > 0.05eV^. Classical /3-decay experiments pro- 
vide an upper limit on the mass of the electron neutrino of 



fixing the Hubble constant implies ^15 eV (Gershtein 



2.2 eV (Lobashev 2003 1, but may in future reach constraints 
as low as 0.2 eV Wolf fc KATRIN collaboration|20"T0 |. 

Sensitivity to well below 1 eV is necessary to distin- 



al. 



20101, 



guish the neutrino mass hierarchy (Jimenez et 
which is of particular interest when formulating theories of 
particle physics beyond the standard model. Fortunately, 
such strong constraints can be obtained from the cosmic 
neutrino background. Simply requiring Q^, ^ JIm ~ 0.3 and 
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fc Zel'Dovich 1966 Lesgourgues fc Pastor 20061. The cos- 



mic microwave background (CMB; e.g. Komatsu & et al. 
(20111) combined with measurements of the matter power 



spectrum from galaxy surveys provide more stringent con- 



straints, presently around M„ < 0.6 eV (Elgaroy et al. 


2002 


Seljak et al. 


2006 


Tegmark et al.|2006||Gratton et al. 


2008 


Mantz et al. 


2010 


Thomas et al.|2010 Reid et al.|2010 


1. Sel- 


jak et al. { 2006 1 combined the Lyman-a forest flux power 



spectrum with the CMB to claim the tightest constraint 
so far; Af^ < 0.17 eV. This strong limit may, however, be 
partly caused by a slight tension between the power spec- 
trum amplitude preferred by the Lyman-a forest data and 
that preferred by the CMB. In any case, it should be tested 
and improved by upcoming data ( Slosar fc et al.|2011 1. The 
Lyman-a data alone provides a limit of < 1 eV (Viel 
2005 Viel et aL]|2010 1, slightly stronger than that provided 
by the CMB data alone ([Larson fc et aL||201lK 



The sensitivity of the matter power spectrum to neu- 
trino mass is due to the free-streaming of the neutrinos 
which suppresses the growth of structure. On large scales, 
where matter perturbations can be described by linear the- 



© 2011 RAS 



2 S. Bird et al. 



ory, this effect is relatively straightforward to quantify with 
numerical simulations or approximate schemes (e.g 
(19801) 



et al. 



( 1994 1 ; Valdarnini et al.\{ 1998|) using standard Boltzmann 



Zeldovich et al.| (19801; Ma fc Bertschinger 



Bond 



codes ( Lewis et al.|2000[ Lesgourgues fc Tram|201l" I . Prob- 



ing the effect of neutrinos on smaller, more non-linear scales 
is more challenging, but has been attempted by, for example, 



Saito et al. ( 2008 2009 2011 1, Wong ( 2008 1 and Lesgourgues 



et al. (20091, with extensions to perturbation theory. Also 
of interest here is the approach of [Abazajian et "aL (20051; 
Hannestad et al. (2005), who described the fully non- linear 
regime by perturbatively evaluating the effect of massive 
neutrinos on halo profiles. 

However, accurately probing the effect of massive neu- 
trinos on the growth of structure arbitrarily far into the non- 
linear regime requires fitting formulae precisely calibrated 
from high dynamic range N-body simulations. The formula 



we present here builds on the widely-used HALOFIT ( Smith 
|et al.|2003[ ). HALOFIT models the non-linear matter power 
spectrum by assuming that the growth of structure is due 
to a combination of the density structure of individual DM 
haloes and the clustering of these DM haloes relative to 
each other. It has been calibrated against a wide range of 
CDM simulations, superseding earlier prescriptions based 



on the stable clustering hypothesis (Hamilton et al. 1991 
Peacock & Dodds 19961. At low and moderate redshifts 



HALOFIT models the non-linear growth of structure in a 
ACDM cosmology with an accuracy of < 5 — 10% (Heit- 
mann et al. 20101 for k < l/iMpc~^, but it has not been 



calibrated for non-zero neutrino masses. 

Fortunately, the quality of the numerical simulations 
has advanced to the point where it is possible to model the 
impact of neutrinos to the level of a few percent, despite 
the difficulties presented by the large thermal dispersion of 

' |Viel et al.|[20T0l ). A 



light neutrinos (Brandbyge et al.^^2008 



quantitative numerical study of the suppression of the mat- 
ter power spectrum due to free-streaming of neutrinos ap- 
pears timely, as future matter power spectrum surveys aim 
to make measurements accurate to a few percent. Realising 
the potential of these measurements to constrain the neu- 
trino mass requires tight control of systematic error from 
theoretical predictions. 

We present results from simulations with the re- 
cently developed version of the Tree-SPH code GADGET -3 
(Springel 2005 Viel et al. 20101. This suite of simulations 



was designed to evaluate the performance of HALOFIT pre- 
dictions when compared to full numerical simulations. We 
use our simulations to calibrate an improved fitting formula, 
capable of accurately reproducing the effect of massive neu- 
trinos on the non-linear matter power spectrum. We discuss 
our methods in Section [21 present our results in Section |3j 
and conclude in Section HI 



2 MODELLING THE NON-LINEAR MATTER 
POWER SPECTRUM 

Here we briefly describe the numerical methods we use to 
simulate structure growth in the presence of massive neutri- 
nos, as well as the approximations underlying HALOFIT . In 
Section [2.1| we review the predicted effect of neutrinos in lin- 
ear theory. In Section [2. 2| we discuss the strategies used for 



Matter Power Spectra at 2:^0,1 




10" 
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Figure 1. Matter power spectra from a simulation with mas- 
sive neutrinos (dashed lines), and A/^ = 0.6 eV, compared to 
the corresponding power spectra from a ACDM simulation (solid 
lines). Dotted lines show the linear theory prediction for mas- 
sive neutrinos, while dot-dashed lines show the linear theory pre- 
diction for ACDM. Upper (green) lines show spectra at z = 0, 
while lower (black) lines show them at 2 = 1. The simulations 
shown are S60 (small scales) and L60 (large scales), whose pa- 
rameters are described in Table ^ S60 shows only modes where 
k > 4 X 27r/ Box-size, so sample variance is small, and L60 is 
cut off at k = OAhMpc~^ , for clarity. Note that the absolute 
power spectrum is slightly less converged with respect to box size 
than the ratio of the massive neutrino power spectrum with and 
without massive neutrinos. 



simulating the effect of massive neutrinos, describing meth- 
ods based on both particle and fourier-space representations 
of neutrinos. In Section |2.4| we summarise HALOFIT . Fi- 
nally, Section [2 . 3| gives details of our numerical simulations. 



2.1 Massive neutrinos in linear theory 

Neutrinos in the very early Universe are relativistic and 
tightly coupled to the primordial plasma. They decouple 
while still relativistic and then redshift adiabatically. The 
transition to non-relativistic behaviour takes place when 
PL = 3n, at a redshift of {1 + z) = 2 x 10^(M^/1 eV). How- 
ever, it is some time later before relativistic corrections are 
negligible, especially in the high-velocity tail of the Fermi- 
Dirac distribution. Once non-relativistic, massive, weakly- 
interacting neutrinos behave qualitatively as a species of 
warm/hot dark matter, suppressing fluctuations on scales 
smaller than their thermal free-streaming length. 

This suppression has a distinctive shape in the linear 
regime. Neutrinos of mass Af^ have energy density 57^, and 
make up a fraction, of the dark matter, where 



93.14 /i2eV ' 



(1) 



Qm is the present day matter density and h is the Hubble 
constant. The residual comoving peculiar thermal velocity 
of the neutrinos at redshift z is 



150(1 



leV 



km/s . 



(2) 



By analogy with the Jeans length, we can deflne a co-moving 
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free-streaming scale 



(3) 



2 H{t)a{t) 

3 vt\i(t) 

Neutrinos cannot cluster on scales smaller than fcps, and so 
the matter power spectrum, defined as 

P{k) = \~5{k)\\ (4) 

where a tilde denotes a Fourier transformed quantity, is sup- 
pressed by a constant factor proportional to the neutrino 
energy density, fv The Boltzmann equation can be solved 
numerically to show that the fractional change in the power 
spectrum for < 0.07 is 

On larger scales, k < kps, neutrinos have free-streamed only 
for a portion of the cosmic history, inducing less suppression. 
The largest scale on which neutrinos alter the matter power 
spectrum is given by the wavenumber for which they free- 
streamed only immediately after becoming non-relativistic 
I Lesgourgues fc Pastor||2006[ ) 



fcnr ~ O.OlSWfiM 



leV 



/iMpc"^ , 



(6) 



Today, neutrinos in the observationally acceptable range 
have fcnr ~ 10"'^, placing the free-streaming effect at the 
scales probed by galaxy and Lyman-a measurements of the 
matter power spectrum, as illustrated by Figure [l] which 
shows non-linear matter power spectra estimated from our 
simulations. For further information on the effect of neutri- 
nos in linear theory see, e.g. Kolb & Turner (19901; Bond 



et al. ( 19801; Eisenstein & Hu 



( 1999 1 and references therein. 



2.2 Incorporating massive neutrinos into N-body 
simulations 

The first N-body simulations incorporating massive neu- 
trinos were aimed at evaluating their potential as a dark 



matter candidate (White et al. 1983 Klypin et al. 1993 



|Ma fc Bertschinger 



1994[ ). Recently attention has shifted to 



their perturbative effect as a sub-dominant component in a 
mainly ACDM universe. We build on the work of |Viel et al.| 
( |2010p, who inco rporated massive neutrinos into GADGET - 
3 ( |Springel|2005 l in two ways; as a new particle species, and 



using a Fourier-space representation ( Brandbyge & Hannes- 
|tad|2009[ ). |Viel et alT] ( |2010[ ) have described the implementa- 
tion in detail and performed extensive validation tests. Here 
we briefly review the two implementations of neutrinos. We 
refer the reader to [Brandbyge et al.] (|2008[ ) ; [Brandbyge fc| 
Hannestad[ ((20091 and[Viel et al.[([2010|) for further details. 



Both the particle and the Fourier space implementations 
incorporate massive neutrinos directly into the simulations 
and so are significantly more accurate than simulations al- 
tering only the initial conditions of a ACDM simulation, as 
e.g. in Agarwal fc Feldman (20111. 

All our simulations are inherently Newtonian, ignoring 
relativistic corrections. It is therefore important to set up 
the initial conditions of the simulation at sufficiently low red- 
shift that the vast majority of neutrinos are non-relativistic. 
Otherwise, neutrinos will be treated as massive while still ef- 
fectively massless, over-estimating their effect by a constant 
factor in the linear regime, at the highest redshifts. 



2.2.1 Particle-based Neutrinos 

Massive neutrinos can be incorporated into cold dark mat- 
ter simulations as a separate low-mass collisionless particle 
species. Structure growth is primarily driven by the cold 
dark matter. The neutrino species has significant thermal 
dispersion, preventing it from clustering on small scales. 
This allows us to avoid following the small-scale evolu- 
tion of the neutrino component, first by ignoring the short- 
range tree force for the neutrinos and second by setting the 
timestep purely from the dark matter species, effectively re- 
laxing the Courant conditions for the neutrino component. 



Brandbyge et al. ( 2008 1 and Viel et al. (20101 found that 



these approximations significantly speed up computation 
while having a negligible impact on results. 

One potential problem with particle based simulations 
is thermal shot noise from the finite sampling of the neu- 
trino distribution (Wan g fc White[[2007 1, especially for the 
higher thermal velocities encountered at high redshift or for 
particularly low mass neutrinos. We deal with this by ex- 
plicitly checking the effect of shot noise, and increasing the 
neutrino particle number should it be non-negligible. 



2.2.2 Fourier-space Neutrinos 

The effect of neutrinos on the dark matter can also be ap- 
proximated by adding their gravitational force to that of 



the dark matter in Fourier space (Brandbyge & Hannestad 
12009'), as part of the particle mesh computation. The massive 
neutrino component is assumed to have the power spectrum 
predicted by linear theor y, in our case com puted using the 
Boltzmann code CAMI^ ( Lewis et al.|2000 |. Since the neu- 
trinos are not followed directly, we can avoid shot noise and 
the need to deal with large thermal velocities when com- 
puting their effect on the overall growth of structure. How- 
ever, because the neutrinos are evolved using linear theory, 
Fourier-space simulations neglect any non-linear evolution in 
the neutrino component itself, including back-reaction from 
non-linearities in the dark matter. As we will discuss in Sec- 
tion [Sj these effects are important, and so we shall focus 
mainly on simulations involving particles. 



2.3 The numerical simulations 

We primarily used particle-based simulations, summarised 
in Table[l] We also ran several simulations using the Fourier- 
space implementation of neutrinos using a grid with 512"^ 
elements, but as these were not used in the final results we 
omit them from the table for clarity. 

Our fiducial simulations have 512^ dark matter parti- 
cles in a box of 150Mpc/i~\ We tested = 0.15 eV, 
0.3 eV and 0.6 eV, running 6 simulations for each value of 
A'li, with different realisations of the initial conditions, to 
suppress cosmic variance. To check we had correctly sam- 
pled the large scale modes, to explore the linear regime and 



to compare directly to Viel et al. ( 2010 1, we ran also simula- 



tions with box sizes of 512 Mpc h .To facilitate comparison 
withlBrandbyge et al.[( [2008[ ) and [Viel et"ar| ( [2010[ ), our fidu- 
cial cosmology had S7m = 0.3, a spectral index Us = 1.0, the 
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NcimG 


M„ (eV) 


Box (Mpch^^) 


JVl/3 

"dm 


JVV3 

bar 




No. Seeds 




^ M 


As/lO"" 




h 


cTg (derived) 


S15 


0.15 


150 


512 





512 


6 


24 


0.3 


2.27 


1.0 


0.7 


0.84 


S30 


0.3 


150 


512 





512 


6 


49 


0.3 


2.27 


1.0 


0.7 


0.80 


S60 


0.6 


150 


512 





512 


6 


99 


0.3 


2.27 


1.0 


0.7 


0.73 


L15 


0.15 


512 


512 





512 


1 


24 


0.3 


2.27 


1.0 


0.7 


0.84 


L30 


0.3 


512 


512 





512 


1 


49 


0.3 


2.27 


1.0 


0.7 


0.80 


L60 


0.6 


512 


512 





512 


1 


99 


0.3 


2.27 


1.0 


0.7 


0.73 


S15NU 


0.15 


150 


512 





1024 


1 


24 


0.3 


2.27 


1.0 


0.7 


0.84 


S60NU 


0.6 


150 


512 





1024 


1 


99 


0.3 


2.27 


1.0 


0.7 


0.73 


S60AS 


0.6 


150 


512 





512 


1 


99 


0.3 


2.0 


1.0 


0.7 


0.73 


S60NS 


0.6 


150 


512 





512 


1 


99 


0.3 


2.27 


0.9 


0.7 


0.73 


S60ND 


0.6 


150 


640 





640 


1 


99 


0.3 


2.27 


1.0 


0.7 


0.73 


S60H 


0.6 


150 


512 





512 


1 


99 


0.3 


2.27 


1.0 


0.75 


0.79 


S30OM 


0.3 


150 


512 





512 


2 


49 


0.25 


2.27 


1.0 


0.7 


0.68 


S30OL 


0.3 


150 


512 





512 


1 


49 


0.4 


2.27 


1.0 


0.7 


1.07 


S60BA 


0.6 


150 


512 


512 


512 


1 


99 


0.3 


2.27 


1.0 


0.7 


0.73 



Table 1. Summary of simulation parameters, erg has been derived from the other parameters. For each simulation with massive neutrinos 
we ran a ACDM simulation with identical parameters, but without massive neutrinos, in order to compute their ratio. The value of erg 
given is that for the simulation with massive neutrinos. 



Hubble parameter h — 0.7, and amplitude of the primor- 
dial power spectrum Ab = 2.27 x 10~^. This implies that a 
simulation without neutrinos has as = 0.87 at z = 0, but a 
simulation with A/„ — 0.3 eV has ag — 0.8. We tested our 
formula using a variety of simulations with different cosmo- 
logical parameters, as summarised in Table [l] 

The initial redshifts of our simulations were Zi — 99 for 
= 0.6 eV, Zi = 49 for = 0.3 eV and z^ = 24 for 
Ml, = 0.15 eV, which we found to be the earliest redshift 
where the neutrinos were fully non-relativistic. Our initial 
conditions were generated using the Zel'dovich approxima- 
tion, and the linear theory power spectrum was calculated 
by CAMB. We acknowledge that second order Lagrangian 
perturbation theory is in principle more accurate, but our 
simulations are started at sufficiently high redshift that the 
effect should be small. The transfer function used for our 
cold dark matter (CDM) particles was a weighted average 
of the linear theory transfer functions for CDM and baryons, 
to account for the slight difference between them ( |Somogyi| 
& Smith 20101. Neutrinos and CDM had identical initial 



random phase information, to simulate adiabatic initial con- 
ditions. 

The numerical convergence of our results and the im- 
pact of neutrino shot noise were checked with simulations 
S60ND, S15NU and S60NU. We further checked the effect 
of baryonic physics with simulation S60BA, which included 
512^ baryon particles. S60BA modelled radiative cooling and 
star formation using the default multiphase model of GAD- 
GET -3 ( [Springel fc Hernquist|2003| . Differences much larger 
than 10% have been found for fe > 20 /iMpc~^ between simu- 
lations implementing different radiative processes (e.g. metal 



cooling) or feedback recipes in ACDM models (see e.g. Rudd 
et al.| ( |2008| ); [Guillet etl[L| ( |2010[ ); [van Daalen et al.| ( |2011[ )) 



Furthermore, the model of AGN feedback suggested by |van 



Daalen et al. (20111 in order to reproduce observations of 



groups of galaxies produces a power spectrum suppression of 
roughly 10% at fc = 1/iMpc"^ and 1% at fc = 0.3/iMpc"^ 
We did not include this effect; however, [van Daalen et aL] 
(20111 found it to be independent of cosmological parame- 



ters, suggesting that it is also independent of the neutrino 
mass, and thus should not affect our results. 

2.4 The HALOFIT formula 



HALOFIT ( Smith et al. 2003 1 estimates non-linear growth 
from a given linear theory power spectrum, based on the 
principles of the halo model. These assert that the growth 
of halos should depend only on the local physics at the scale 
of the halo, and not the details of the pre-coUapse mate- 
rial, nor the larger-scale distribution of matter. Thus growth 
deep in the non-linear regime becomes independent of the 
details of the linear power spectrum, depending only on the 
non-linear scale and the slope and curvature of the power 
spectrum. In the mildly non-linear regime, the non-linear 
correction is assumed to be a function of the instantaneous 
linear power spectrum, without explicit dependence on the 
historical growth rate. HALOFIT is accurate to 5 — 10% at 
z < 3 ( [Heitmann et al.|2010[ ), suggesting that these approx- 
imations are reasonably accurate for cold dark matter. 

HALOFIT realises the above by splitting the non-linear 
power spectrum into two components. A quasilinear term 
Aq, and a non-linear term Ah. This means the dimension- 
less power spectrum can be written as, 

^ ' =Al + Al. (7) 



27r2 



+ A|. 



Both terms depend on the non-linear scale, k^, defined as 
the wave-number where the perturbation amplitude, filtered 
by a Gaussian, becomes greater than unity and linear theory 
breaks down. 

The quasilinear term, Aq, is given as an enhancement 
to the linear power spectrum, and incorporates an exponen- 
tial cut-off ensuring it is negligible on small scales. The form 
of the quasilinear term is, 



A|(fc) = AE(fc) 



(1 + A?,(fc))'3(") 



l + a{n)Al{k) 



Here 3 + n — 



dlna^ 
dlni? 



exp[-/(y)]. (8) 



(9) 
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so that n is the effective spectral index of the power spec- 
trum. /3 and a are functions of n and f{y) — y/4-|- 3/^/8. We 
refer readers inter ested in the exact v alues of these functions 
to Appendix C of Smith et al. (20031. is the dimension- 



less linear theory power spectrum. 

Non-linear growth, Ah, dominates when y = k/ka- 
The form of Ah used by HALOFIT is. 



A|(fc) 



A|'(fc) 



Where C 



1 + li{n)/y + u(n)/y^' 

a{n,C)y^^''-'''> 



> 1. 



(10) 



l-h6(n,C)y-f2(") -I- [c(rz,C)/3(J7)y]a-T(".c) ' 

(11) 

if In a^{R) 



(12) 



the spectral curvature of the power spectrum, a, b, c, fj, and u 
are functions of n and C, and vin) is not related to neutrinos. 
f 1,2,3 are functions of Qm, the total matter density. For the 
values of these functions, we again refer to Appendix C of 



Smith et al. ( 2003 1. Note that the non-linear growth depends 



only on y, n, C and the matter fraction. 

Massive neutrinos affect this prescription because they 
suppress the linear power spectrum, increasing the non- 
linear scale and retarding growth. HALOFIT estimates the 
effect of this suppression on non-linear power, but neglects 
any effect of the neutrinos on the process of halo collapse, 
as well as the effects of non-linear growth on the neutrinos 
themselves. 



3 RESULTS 

Figure [2] shows ratios between the matter power spectrum 
for massive neutrinos, and ACDM. Plotted are the predic- 
tions of linear theory, of simulations S60 and L60, and of 
grid-based simulations with identical parameters. The max- 
imal power spectrum suppression in linear theory is propor- 
tional to —8fv, compared to a maximal non-linear suppres- 
sion of approximately —10/^ for z ^ 1, and —10.5/,/ for 
z = 2. As discussed in Section [2. 1[ linear theory predicts an 
increasing suppression approaching an almost flat plateau on 
small scales. In the numerical simulations, the flat plateau 
is replaced by a minimum. The suppression increases up 
to ~ l/iMpc~^, after which it gradually decreases. This 
distinctive spoon shaped pre-virialization feature was also 



found by Brandbyge et al. 



and in CDM simulations by 



(20071 



(20081 and 



Viel et al. (20101, 



Lokas et al. ( 1996 1; Smith et al 



We can understand this shape by considering the effect 
of neutrinos on non-linear growth. The suppression caused 
by massive neutrinos delays the onset of non-linear growth 
and increases the wave-number of the non-linear scale. Be- 
tween the non-linear scale for the simulation without neutri- 
nos and the non-linear scale for the simulation incorporating 
neutrinos, the suppression will be enhanced by the absence 
of non-linear growth. On the smallest scales, the non-linear 
growth lost as a consequence of the massive neutrinos be- 
comes an ever-smaller fraction of the total, and the extra 
suppression decreases. Eventually non-linear growth comes 
to dominate the linear effect, and the suppression is less than 
predicted by linear theory. 



As mentioned in Section |2.3| we performed checks for 
numerical convergence. A simulation with increased reso- 
lution, S60ND (see Table [T]), had a relative matter power 
spectrum which differed from the fiducal simulation, S60, 
by 1 — 2% for k < 7/iMpc~^. Furthermore, the relative mat- 
ter power spectra for our two different box sizes, 512 and 
150Mpc/i~^, were always in good agreement on the largest 
scales probed by the smaller box, showing that our smaller 
box size was sufficient for an accurate estimate of the sup- 
pression of matter power spectrum due to the free-streaming 
of the neutrinos. The larger box showed a slightly increased 
suppression on small scales, probably due to limited reso- 
lution. We found that the effect of baryons was, at all red- 
shifts, less than 1% at < 8/iMpc~^, gradually increasing 
on smaller scales, in agreement with the results of e.g. |Jing| 
al.| ( [2006l ); |van Daalen et al.| ( [20TT| . 
Simulations S15NU and S60NU were run with a higher 
number of neutrino particles, to check explicitly the effect 
of shot noise. We found that, although the neutrino power 
spectrum is indeed shot noise dominated on small scales, 
the power on these scales is so small that their impact on 
the growth of the dark matter power spectrum is negligible. 
Hence shot noise has only a minor impact for simulations 
with a box size of 150 Mpc at the level of 1% at the 
relevant scales and redshifts. 



3.1 Particle and Fourier-space implementation of 
neutrinos 

Figure [2] shows that particle and Fourier-space implementa- 
tions of neutrinos give similar results for the matter power 
spectrum, and both agree with linear theory on large scales. 
The particle implementation, however, shows a smaller sup- 
pression, which scales roughly linearly with M,^, particularly 
near the trough of the dip. 

This discrepancy has again a simple physical explana- 
tion; structure growth in the dark matter induces structure 
growth in the neutrino component, as the less energetic neu- 
trinos fall into the gravitational wells created by the dark 
matter. This increases the power in the neutrinos, drawing 
them closer to the non-linear dark matter power spectrum 
and leading to a reduced neutrino suppression. Essentially 
this is a backreaction effect, where the dark matter drags 
the neutrinos with it; since it results from the non-linear 
growth, it is not being fully accounted for by linear theory 
neutrinos. Note that structure growth in the neutrino com- 
ponent itself remains linear; the effect is due to non-linear 
growth in the dark matter. Because this effect is important 
at the accuracy we wish to achieve, our main results (Fig- 
ures |3j [4] and |5| are based on simulations with the particle 
implementation of neutrinos. 



3.2 Comparison to the HALOFIT model 

Figures [3] |4] and [5] show the main results of this paper. 
Each Figure shows the suppression of the matter power spec- 
trum caused by massive neutrinos at four redshift snapshots 
{z — 0, 0.5, 1 and 2), compared to the predicted effect from 
HALOFIT and linear theory. Figure [3] shows = 0.6 eV, 
Figure |4] has = 0.3 eV and finally Figure [5] has = 
0.15 eV. 
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Figure 2. The ratio between the matter power spectrum at z = (Left) and 2 = 1 (Right) for simulations with and without massive 
neutrinos {Mi, = 0.6 eV). Solid green lines show results for simulations using particles, while dashed red lines show results from Fourier- 
space simulations. The black dot-dashed lines show the predicted effect from linear theory. The simulations used were L60 and S60; box 
sizes were 512Mpch~^ and 150Mpc/i~^, with parameters shown in Table ^ For the smaller box sizes, more than one realisation of 
structure is available, and hence we show (dotted grey lines) the one a error due to sample variance, as estimated numerically. Particle 
and Fourier-space methods use the same initial structure realisations. The Fourier-space method has more power in the largest scale 
mode, probably because the phase structure of the neutrino component has been prevented from evolving, enhancing the effect of sample 
variance. 
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Figure 3. The effect of massive neutrinos on the matter power spectrum for a neutrino mass of Af^ = 0.6 eV. Solid lines show the 
ratio between simulations with and without massive neutrinos, for both L60 (red), with a 512Mpc/i~^ box and S60 (orange), with a 
150Mpc?i~^ box. Initial redshift was 99. The blue dashed line shows the estimated ratio using HALOFIT , while the black dashed line 
shows the prediction from linear theory. 
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Figure 4. The effect of massive neutrinos on the matter power spectrum for a neutrino mass of Mi, = 0.3 eV. SoUd lines show the 
ratio between simulations with and without massive neutrinos, for both L30 (red), with a 512Mpc/i~^ box and S30 (orange), with a 
150Mpc/i~^ box. Initial redshift was 49. The blue dashed line shows the estimated ratio using HALOFIT , while the black dashed line 
shows the prediction from linear theory. 



HALOFIT clearly over-predicts the suppression of the 
matter power spectrum due to massive neutrinos in the non- 
linear regime. The cause is similar to that discussed in Sec- 
tion [STT] HALOFIT includes massive neutrinos only through 
the linear theory neutrinos suppression on the non-linear 
scale, and thus neglects any back-reaction from the dark 
matter. This interpretation is supported by the good agree- 
ment of HALOFIT with the Fourier-space simulations. 

The largest discrepancy, around 10% of the total sup- 
pression, occurs at ~ 1 ftMpc"^. The location of the maxi- 
mal suppression in the numerical simulation moves to larger 
scales at lower redshifts, an effect which is again not cap- 
tured by HALOFIT , although it was present to some extent 
in our Fourier-space simulations. Furthermore, the ampli- 
tude of the suppression decreases with redshift, which is to 
be expected from our discussion in |3.1[ as non- linear growth 
occurs in the dark matter more neutrinos fall into the grav- 
ity wells. A cross-over redshift occurs at 2 = 1; on very 
small scales the suppression here is equal to that of linear 
theory, while at lower redshifts it is less. Note also that in 
the quasilinear regime, 0.05 < k < 0.2/iMpc~'^, the simula- 
tions clearly agree much better with HALOFIT rather than 
linear theory The dependence of our results on Mi, is well 
described by a linear relation, with the maximal suppression 



for a given redshift being proportional to f,^, although the 
redshift dependence is more complicated. 

HALOFIT over-predicts the effect of neutrinos on the 
smallest scales. This is not due to neutrino physics, but is 
a discrepancy induced because, as found by |Hilbert et al.| 
( 2009 1 , HALOFIT under-predicts the growth of non- linear 
power for fc > 2 /iMpc"^ in a ACDM universe by up to a fac- 
tor of two. We corrected this by re-fitting the HALOFIT pa- 
rameter that controls the small-scale power using our ACDM 
simulations, after which we could reproduce the asymptotic 
neutrino effect. We detail our modifications in Appendix [A] 
Note that HALOFIT also has reduced accuracy at z = 3 — 4, 
with errors of 15 — 20%, compared to 5 — 10% at z; = — 2. 
Because of this, it fails to accurately predict the location of 
the peak non-linear suppression at z > 2. We did not correct 
this effect, as our attempts were found to negatively impact 
accuracy at lower redshifts. 

We have also performed some simulations varying the 
cosmological parameters from our fiducial values, as de- 
scribed in Section [22] The results of these simulations were 
similar to those for our fiducial cosmology, and agree with 
the results discussed in Viel et al. ( 2010) . The dependence of 
the non-linear suppression on cosmology is largely captured 
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Figure 5. The effect of massive neutrinos on the matter power spectrum for a neutrino mass of M,j = 0.15 eV. Sohd hues show the 
ratio between simulations with and without massive neutrinos, for both L15 (red), with a 512Mpc/i~^ box and S15 (orange), with a 
150Mpc/i~^ box. Initial redshift was 24. The blue dashed hne shows the estimated ratio using HALOFIT , while the black dashed line 
shows the prediction from linear theory. 



by HALOFIT , with the exception of a weak dependence on 
Q,M, which we include in our fitting formula. 



3.3 Towards an improved fitting formula 

In the previous Section, we compared our simulations to 
HALOFIT . While there was some discrepancy, it seemed to 
capture many of the important physical features. In hght of 
this, we have designed an improved fitting formula by modi- 
fying HALOFIT to explicitly account for the effect of massive 
neutrinos. Our fitting formula was derived by performing a 
series of least-squares minimisations of £ = Asi,„/A|t — 1, 
where Agj^ is the matter power spectrum from the simula- 
tion and A|t is the equivalent from the fitting formula. Both 
spectra have been rebinned to be evenly spaced in log k. To 
find a minimum, we used Powell's method as implemented in 
SciPjj^ To avoid any contamination from sample variance, 
we omitted scales with k < 6j|^, where Box is the box size. 
The redshift bins used in our fits were 2 = 0, 0.2, 0.5, 1, 2, 3. 
HALOFIT models the effect of non-linear growth at z = 4 
rather poorly, and there is some residual shot noise present 
at these redshifts in our neutrino simulations. Furthermore, 



http : / /www . scipy . org 



the only current probe of the matter power spectrum at 
these redshifts is the Lyman-a forest, which should be ex- 
amined by modelling the flux power directly with full hydro 
simulations, a; 



Viel et al. 



(120101. Hence we do not at- 
tempt to provide a good fit for z > 3. We checked that our 
fit was not significantly altered by excluding the simulation 
results at z = 3. 

Before addressing the effect of neutrinos, we attempted 
to correct the lack of power shown by HALOFIT on the small- 
est scales. Eq. |12| shows that the asymptotic behaviour of the 
non-linear term in HALOFIT is given by 



a{n) 



3(/i(n)-l)+7(n) 



y 



,7(") 



(13) 



Thus, to correct the behaviour on small scales we altered 
j{n), fitting it to the ACDM simulations labelled in Table[l] 
as S60, S60AS, S60NS, S60H and S30OL. We were interested 
only in the asymptotic behaviour, so we omitted the results 
from the large simulation boxes from the fit. We thereby 
included extremely small scales into the fit; down to k — 
30hMpc~^. The numerical minimisation algorithm proved 
unstable when constrained to < 6/iMpc~'^. We checked 
that the resulting fit is a good match to larger scales. 

The largest remaining contribution to the discrepancy 
between HALOFIT and our simulations is the decrease in the 
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maximal neutrino suppression due to neutrinos falling into 
dark matter halos. This is proportional to and roughly 
constant on large scales. To account for it we modified the 
nonlinear term with the ansatz 

If" 



1 + mk^ ' 

Anl(1 + Q.): 



(14) 
(15) 



with / and m found by minimisation. We considered other 
powers of k in the denominator of Qv, as well as additional 
terms and powers of the non-linear scale ka, but none of 
these significantly enhanced the fit. 

To correctly match the location of the maximal sup- 
pression at higher redshift, z ^ 2, we needed to alter the 
quasilinear term, refiecting a leakage of non-linear power 
from small scales to larger scales caused by the high ther- 
mal dispersion of the neutrinos. Therefore, we modified Eq. 
[9] as follows, 

' {1 + Al{k)f^"^' 



Ai 



1 + 



l + a{n)Al{k) 
1 + 1.5fe2 



/3(n) ^ P(n) + f,{r + sn^) , 



exp[-/(j/)], (16) 

(17) 
(18) 

where p, r and s were found again by minimisation. These 
coefficients, as well as I and m were found by simultane- 
ously fitting to the simulations S15, S30, S60, L15, L30 
and L60. For the larger box sizes, we used only scales 
with k < 0.5 /iMpc~^, and for the smaller boxes, we used 
k < 6/iMpc~\ 

Finally, there was some residual dependence on Q,m- To 
account for this, we modified Z as / — ^ Z — t(r2jv/ — 0.3), and fit 
t separately, based on a minimisation with simulations S30, 
S30OM and S30OL. The numerical values of the coefficients 
are shown in Appendix [X] together with a summary of the 
fitting formula. The reduced error in the fit can be defined 



£ 



(19) 



where A*'d is the number of data points used in the fit. We 
find £ = 0.007 with our coefficients, compared to £ = 0.04 
with the unmodifed HALOFIT . 

A rough estimate of the error on our fitting formula, 
including realisation noise at the 1 — a level, is given by 

e = E{k,z)P{k), (20) 

- l + logll + fc/fc^-)]-^" ' ^ ' 

where kcr{z) is again the non-linear scale. Differences be- 
tween the fitting formula and the simulations change sign 
between z = and z = 0.2, suggesting the dominant contri- 
bution is realisation noise. The quoted errors are valid for 
k < 7/iMpc~^, and 2 ^ 3. At z > 3 our fitting formula is 
somewhat less accurate. The maximal error over all scales is 
E 

max — /„; for = 0.3, E 

max 

2%. 

Although they were not included in the fit, we checked 
that our formula accurately reproduced simulations with 
different cosmologies, chosen to be at the bounds of the 
range allowed by current data. The fitting formula errors for 
S60AS, S60NS, and S60H were very similar to those found 



when compared to the fiducial cosmology, although there 
were one or two outliers caused by realisation noise. Check- 
ing against S30OM, which has Q,m = 0.25, produced some- 
what larger errors at the highest redshift bin, with z = ?>. 
The lower value of Q.m means that z = 3 is deeper into the 
quasilinear regime, where HALOFIT is less accurate. 

Figure [6] shows the performance of our improved fitting 
formula, in a similar format to that of Figure |3] It is clear 
that our improved formula captures the suppression due to 
neutrinos significantly more accurately than HALOFIT , cor- 
rectly predicting the location and depth of the maximal sup- 
pression to within the quoted errors. 



4 DISCUSSION AND CONCLUSIONS 

We have performed a suite of detailed simulations of the 
matter power spectrum, incorporating the effects of mas- 
sive neutrinos using both particle and Fourier-space meth- 
ods, but focussing on simulations with neutrino particles. 
We have compared our results to those predicted by the 
widely used HALOFIT approximation to the non-linear mat- 
ter power spectrum, and presented an improved fitting for- 
mula. 

Observational constraints on the neutrino mass prior to 
this work have either used linear theory ( Reid et al.|2010 1 or 
HALOFIT to estimate the effect of massive neutrinos on the 



matter power spectrum. Thomas et al. (20101, for example, 
used photometric redshifts at z = 0.45 — 0.65, on scales 
k < 0.2/iMpc~\ This is in the weakly non-linear regime, 
where we find HALOFIT to be accurate. Hence our results 
should not significantly change their published constraints 
and that of studies probing similar scales. 

However, while HALOFIT reproduces the suppresion of 
the matter power spectrum due to the free-streaming of neu- 
trinos well on some scales, in the fully non-linear regime it 
systematically over-predicts the level of the effect. Further- 
more, HALOFIT fails to capture the redshift dependence of 
the amplitude and scale of the peak suppression. Our quan- 
titative numerical study of the non-linear regime will thus 



certainly be of importance for future galaxy (Wang et al 
Mis': Schlcgol ct al. 2009: C arbone et al.||2011| ), CMB lens- 
ing (^Cooray,^1999, ,Kaplinghat et al.||2003[), Lyman-g forest 



([Gratton et al.||2008[) a nd weak lensing (|Ichiki et aL |2009| 
[Beaulieu fc et al.||2010| ) surveys. [Abazajian fc et al.||201li 
have published a review of the potential constraints on the 
neutrino mass from a variety of observational measurements, 
many of which reach the ~ 0.1 eV range. Realising these 
limits will require theoretical predictions of the power spec- 
trum accurate to the percent level, significantly exceeding 
the accuracy of current prescriptions. 

We have proposed an improved fitting formula which 
modifies HALOFIT to take account of the effects of massive 
neutrinos for k < 7 /iMpc~^ at z ^ 3. Errors are propor- 
tional to the amplitude of the suppression, and thus to 
for A/^ — 0.3 eV, the largest error is roughly 2% of the mat- 
ter power spectrum. The maximal non-linear suppression we 
find is around —10/^, significantly greater than the —8/1/ ex- 
pected from linear theory. Hence constraints obtained using 
our formula should be tighter than predicted by forecasts 
using linear theory. 



Saito et al. (2008 2009 20111 and Wong (20081 have 
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Figure 6. The effect of massive neutrinos on the matter power spectrum for a variety of neutrino masses. Solid lines show the ratio 
between simulations with and without massive neutrinos, for both 512Mpc/x~^ (red) and 150Mpc/i~^ (orange) boxes. The blue dashed 
line shows the estimated ratio using our improved fitting formula, while the black dashed line shows the prediction from linear theory. 



proposed an approximation to the non-linear matter power 
spectrum based on one- loop perturbation theory, while |Les^ 
gourgues et al. ( 2009 1 have modelled the effects of neutrinos 



with an approach based on a time renormalisation group 
flow. These approaches are to some extent complementary to 
our efforts; while these formulae are designed for the quasi- 
linear regime, for k < 0.3/iMpc~^ or z > 1, our formula 
should be useful even in the fully non-linear regime, at the 
lowest redshifts. 

Our fitting formula should be sufficiently accurate to 



make forecasts for the accuracy of neutrino mass measure- 
ments of future surveys. We should stress, however that we 
have concentrated here on the relative impact of the neu- 
trinos; HALOFIT itself only predicts the non-linear power 
spectrum to an accuracy of 5 — 10% i Heitmann et al.|2010 l 
with larger discrepancies at very small scales and high red- 
shifts (2: > 3). Detailed quantitative analysis of actual data 
reaching into the non-linear regime will thus still need care- 
ful calibration with accurate numerical simulations. These 
will have to be tailored to the redshift and scales relevant 
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for the particular data set at hand as the simulations still 
struggle to provide the required dynamical range and the 
fitting formula are not yet accurate at the percent level over 
the full range of scales, cosmological parameters, neutrino 
masses and redshifts of interest. 

We have implemented our fitting formula in the form of 
a series of patches to CAMB, which have been included in 
the publicly released version. 



incorporating the residual dependence on Qm. Finally, the 
altered quasilinear term is 



(1 + A^(fe))«") 



1 + 



l + a{n)Al{k) 



26.3/„fe^ 



exp[-/(j/)] (A8) 



1 + 1.5fc2^ 
/3(n) = /?(n) + /4-6.49 + l.Un^) , 

completing our fitting formula. 



(A9) 
(AlO) 
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APPENDIX A: HALOFIT-;/: AN IMPROVED 
FITTING FORMULA 

The HALOFIT non-linear term is given by 

Ai,'(fc) 



Aiik) = 
Here 3 + n 



(XI) 

a(n,C)j/^-fi'"' 

1 + b(n, C)yMi'^) + [c(n, C)/3(f^)j/]»-^("-C') ' 

(A2) 



din a" 



dlni? 



and C - 



In (7^ 



dlnJ?2 



(A3) 



To allow HALOFIT to more accurately model the growth of 
non-linear power on the smallest scales, we modify it by 
altering 7(72, C) to account for the extra growth on small 
scales as 

7 7 -I- 7' (A4) 
7' = 0.316 - 0.0765n - 0.835C (A5) 

To account for the effect of non-linear growth in the neutrino 
component, the non-linear term is modified as 



/^(2.080 -12.4[»A/-0.3]) 
1 + 1.20 X 10-3fc3 



(a;;l)' = aL(i + Q.): 



(A6) 
(A7) 
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